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Using these models' equivalence to the theory of a scalar field coupled to a linear combination of the 
invariants, we investigate the properties of the propagating modes. Even in the case for which the 
fourth derivative terms in the field equations vanish, we find that the second derivative terms can 
give rise to ghosts, instabilities, and superluminal propagation speeds. We establish the conditions 
which the theories must satisfy in order to avoid these problems in Friedmann backgrounds, and show 

(^ • that the late-time attractor solutions generically exhibit superluminally propagating tensor or scalar 

O ■ modes. 
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I. INTRODUCTION 

The strong observational evidence for an accelerating universe [1-10] has sparked a 
widespread search for a dynamical explanation. Beyond a bare cosmological constant, a plethora 
of other models have been proposed, with quintessence - a dynamical scalar field that behaves 
essentially as a modern day inflaton field, - being perhaps the simplest example (see [11-14]). 
In this context, many potentials have been introduced that yield late-time acceleration and 
tracking behaviour (see [15-24]). 

Among other approaches. Modified Gravity Models have attracted great interest (see [25- 
37]) but also some criticism, partly because they were introduced as purely phenomenological 
models, but more seriously because it was not clear that they possessed a satisfactory Newtonian 
limit in the solar system, or that they were free of ghosts (see [38-43]). 

In this paper, we investigate the propagating degrees of freedom of the so-called CDDETT 
model [26]. There already exist detailed studies of the Newtonian limit [44] and the Supernovae 
contraints [45] for this model. Here we derive conditions that they be free of ghosts, and that 
they have a real propagation speed less than or equal to that of light. 

As we review below, a transformation of the action shows that Modified Gravity models are 
equivalent to a number of scalar fields linearly coupled to higher order curvature invariants. In 
the case in which these curvature invariants are fourth order, the relevant one for the Modified 
Gravity models of Refs. [25, 26], we obtain conditions for the propagating degrees of freedom to 
be well-behaved in their late-time attractor solutions (Friedmann-Robertson Walker spacetimes 
with accelerating expansion). This extends previous work which established their consistency 
in de Sitter backgrounds [38, 39, 44]. 

We find that while untroubled by ghosts, the accelerating power-law attractors in general 
have superluminal tensor and scalar modes, which may place severe theoretical constraints on 
these models. 

II. THE PHYSICAL DEGREES OF FREEDOM 

Our starting point is the action proposed in [26], which we write in the form 

S = j d'x,r9YlR-e^li'-+^f{Z,)], (1) 

where 7 is a constant, Zi = i?^, Z2 = Rai3 R"^ and Z^ = Rap^s R^^"'^ . We have introduced 
9^ = sign(yu) for generality, but note that its presence does not change the late time behaviour 



(3) 



of the accelerating attractors, since for an accelerating universe both the R (Einstein-Hilbert) 
term and the dark matter density become negligible (in other words the exponent of the power 
law attractor does not depend on /i, see [26]). Finally, we take the function f{Zi) to be of the 
form 

^(^•' - (^ • '^) 

where a sum over i is implied. 

The action (1) can be written as that of Einstein gravity coupled to a scalar field, a form 
more suitable for analysing the propagating degrees of freedom (see the Appendix for a general 
analysis). Consider 

where, of course, 7^ 0, otherwise the action is not finite. The variation of this action with 
respect to leads to 

(p = aiZi , (4) 

and, using this relation, action (3) and action (1) yield the same equations of motion. Note 
that when 02 = —4 03 and ai = a^, this action is equivalent to Einstein-Hilbert gravity coupled 
to a single scalar through a Gauss-Bonnet (GB) term Rq^ = R^ — AR^yR'^" + R^yp^R^^f" . 

The coupling of a scalar field with a quadratic expression of the curvature invariants emerges 
naturally in the context of string theory. In particular, as was shown in [46] by Gross and Sloan, 
in the low-energy effective action the dilaton is coupled to a Gauss-Bonnet term. It is well known 
that such a term, expanded about a Minkowski vacuum, ensures that the theory is ghost-free 
(see [47]). 

It might then seem that taking the Oj to be the GB combination is a sensible choice, be- 
cause string theory predicts such a coupling to exist and string theory does not have ghosts. 
However, in models like ours, for which Minkowski spacetime is not a solution, choosing the 
GB combination of parameters Oj is not a sufficient condition for the non-existence of ghosts. 

A ghost is a propagating degree of freedom whose propagator has the wrong sign, and which 
therefore gives rise to a negative norm state on quantisation. Such states are allowed off-shell 
in gauge field perturbation theory, but are unacceptable as physical particles. A theory of 
gravity with fourth order derivatives in the kinetic term inevitably has ghosts [48, 49], but 
even a theory with second order derivatives alone has other potential problems. Once we break 
Lorentz invariance, as in a Friedmann-Robertson- Walker (FRW) background, the kinetic terms 



of a field, even though second order in derivatives, may still have the wrong sign, or may give 
rise to a propagation speed which is greater than 1, or imaginary. To see this in more detail, 
consider the action for a scalar field 0, 

s = J d'xV^ (^^rm' -^sit)v<p'^ . (5) 

The propagation speed of this scalar is v^ = y^S/T. One may wish to impose one or more of 
the following conditions 

1. A real propagation speed: S/T > 0, otherwise all perturbations have exponentially grow- 
ing modes. 

2. A propagation speed less than light: S/T < 1, we will talk about this issue more in detail 
in section III. 

3. No ghosts: T > 0, to ensure a consistent quantum field theory. 

Clearly, unless S and T are positive, and their ratio less than one, we will have instabilities, 
superluminal propagation, or ghosts. We will see that in studying the action for small pertur- 
bations of the metric in Modified Gravity theories we will generally encounter actions of the 
form (5). 

If ai ^ as, the action (1) can be written in terms of an Einstein-Hilbert term plus a particular 
extra piece involving two new scalar fields. Furthermore, because of the special properties of 
the Gauss-Bonnet term, the equations of motion are no longer 4th order, but remain 2nd order 
in the fields. Taking the action (3) and introducing a new scalar field a, we have 

5 = y"rf^Xv/^{[7 + 26a/(0)]i?-f/(0)-6aV(0) + /(0)i?GB} , (6) 

where 

f/(0)=^,/.^"+^^ (7) 



,4n+2 



n 



with b = ai/a^ — 1, and Rq^ the Gauss-Bonnet invariant. 

Making a field redefinition a = x/(26/), the equation of motion for x is x = '^bf'R, and the 
gravitational equations then become 



(7 + X) Rcp - V« V/3X + (7a/3 Dx - 2 i? V^Vpf + 2g^f,Rnf + 8i?(,, V^) V^/ 

,2 



-4R^f,nf -Ag^f.RP'^V.VJ- -AR^^'^^p^V.Vrf -l9ap {l + x) R ~ £j ~ U =0. (9) 



Independent of the background, after application of the Bianchi indentities these equations are 
second-order in derivatives of the fields, thanks to the Gauss-Bonnet combination. 

It is known that adding terms quadratic in the curvature invariants to the Einstein-Hilbert 
action with a cosmological constant yields an extra scalar mode and a spin-2 mode, which is 
generically a ghost because of its fourth-order field equation (see [49-52]). Thus, provided we 
are expanding around a constant background, such as de Sitter space, we can directly infer 
that Modifed Gravity models are generically afflicted by spin-2 ghosts. However, in our case, 
it is clear that the higher derivative terms cancel out identically because 02 = —4 03, as was 
already found in [38, 39, 53]. The remaining extra scalar degree of freedom, x, is due to the 
presence of the extra R^ term in the lagrangian, which vanishes if ai = a^. 

Crucially though, the vanishing of the fourth order term is a necessary but not sufficient 
condition for the absence of ghosts in the spin-2 sector. As we described above, one must also 
separately check the signs of the second order derivatives with respect to both time and space 
(a check was not performed in [53]). 

In this paper we derive and study the kinetic terms for both the spin and spin 2 fields 
in time-dependent backgrounds, to which derivatives of / contribute. We find that both fields 
may be afflicted by instabilities or ghosts, contrary to the claim of the absence of ghosts 
in FRW spacetimes made in [53]. We show that the special case of an empty accelerating 
universe - the late-time attractor of an FRW cosmology in the Modified Gravity model under 
consideration - the propagating states are generically (but not universally) superluminal over 
the (6, n) parameter space. 

III. PROPAGATION IN FRW SPACETIMES 

If the second order derivatives of the spin-2 and spin-0 fields do not have the correct signs in 
FRW spacetimes, then the theory may be inconsistent. The existence of a ghost mode would 
lead, for example, to the over-production of all particles coupled to it. One may think of a 
theory with ghosts as an effective theory, with no ghosts above some cutoff, thereby restoring 
consistency (see [54, 55]). However, this cut-off must be less than about 3 MeV [55]. 

A further condition that one may wish to impose is that the propagation speeds be less 
than or equal to unity. One worry is that the existence of superluminal modes on the rele- 
vant cosmological backgrounds may lead to a catastrophic signature of causality violation (see 
e.g. [56]). Other authors [58-60]) have discussed the problem of superluminal propagation in 



non-lorentzian backgrounds and, in particular, have suggested that the presence of superlumi- 
nal modes would introduce a second horizon, the so called sound horizon, different from the 
light causal horizon, which may lead to ambiguities and inconsistencies in black-hole thermo- 
dynamics [62]. Furthermore, it has been pointed out ([56, 59, 63]) that for some set of initial 
conditions, superluminal modes may yield ill-posed Cauchy problems. In particular, in our case, 
nothing prevents Cq and Cg - the respective speeds of the scalar and tensor modes - becoming 
infinite. In general this would lead to causally connected spatial sections and eventually to an 
ill-posed Cauchy problem. 

On the other hand, in the context of non-commutative geometry, other authors have studied 
superluminal propagation and shown that there is no causality violation if there is a preferred 
reference frame [57]. 

Given these different possibilities, in this paper we will present the constraints from superlu- 
minal propagation in a clearly distinct way from those arising from ghosts, so as to allow readers 
to impose fewer or more constraints, depending on the particular theory they are working with. 

A. Pure Modified Gauss-Bonnet Gravity 

In this section we begin with the special case 

S = j dS^/^ (7i? + f{ml^ - U{<P)) , (10) 

which we refer to as the Modified GB action. We vary with respect to and Qais, write 
^9a/3 = ha/3, SO that Sg""^ = —h""^, and use the many useful identities contained in Ref. [48] to 
derive 



6^S = j d'x^/^g (i?|B/'(0) - U\<P)) 5<P (11) 

5,S = j d'xy^ I ^R [(7 - 4 Of) g-^ + 4V^V"/] 

-i?^^ [(7 - 4 D/) 5^X - A{y,V,f)g^P + K^'V ,f)5f5^^] 
+i?^/(4V^VV) - ^f/(0)^"4 V • (12) 

In order to establish whether the theory is stable and ghost free, we must examine the second 
variation of the action. This can be organized as before, 

5^S = 5lS + 25g5^S + 5lS . (13) 



The easiest term is 

6lS = J d'xV^ {RlJ"{^) - U"{ct>)) 5^ , (14) 

and we can simplify the mixed term using the field equation for 0, to obtain Sg6,j)S = 
J d*x^J—g {6gRQ^)6f, where 6f = f\(j))5(j). After integration by parts we then have 



We simplify the final term using the field equation for gap., and organise according to the number 
of derivatives of hap. In order to check that the original action is ghost free, we need to establish 
that the fourth order terms vanish, and that terms involving two derivatives of the metric have 
the appropriate sign. 

It is already straightforward to see from Eq. (12) that there can be no fourth order derivatives, 
as the terms containing derivatives of the Riemann and Ricci tensors, and the Ricci scalar, have 
already cancelled. The remaining second order terms are 

5,%) = j d'x^^gKA-\[{^ -AUf)g-^ + AV^Vy]{Uh + V'V^h,s) 
+ \[{l- 4n/)5^5,^ - {AV,V,f)g^P + (8V^V^/)5(°(5,^)] 

+ - [4VpV^/] (V^V'/i"^ + V^V^h^f - V^Vr^ - V^Vh'^P) | , (16) 

where A^" = Vuh^"". 

In an FRW space-time, the background fields break Lorentz invariance, with 

g,^u = a^{r)r]p^, (17) 

V.V^/ = 5l5lB{r) + g,.g''C{r), (18) 

where B = f" — 27if', C = Tif (and where we have assumed spatial fiatness). Here a prime 
denotes a derivative with respect to conformal time r, and Ti = a' /a. In this background it is 
convenient to decompose the metric perturbation into scalar, vector and tensor modes in the 
usual way. To check the sign of the kinetic term for the spin-2 particle, we first identify it in 
the expansion of the metric 

hij = gifv + l\ij + C(i\j) + Hij, (19) 



(where the symbol | denotes covariant differentiation with respect to the spatial metric g^- ) as 
the transverse, traceless {H^j\i = 0, H^i = 0) tensor mode Hij. It is straightforward to show 
that the tensor part of the second variation which is second order in derivatives of h^^ is 



^2qT 



j4 / 

a Xy'—g 



^H,, (7 - VO(fi + C)) UW^ + i/J,,(45)V°V°iJ^^^ 



(20) 



Rearranging and using D = VqV^ + A to display the time and space derivatives separately, we 
find 



^aS{2) 



j4 I 

a x^—g 



^H,, (7 - Ag^%B + C)) AW^ + ^H,,{^ - Ag^^C)V^V^W^ 



■ (21) 



We therefore have two conditions for a stable theory free of ghosts: 

7 + 4(/"-W)>0, 7 + 4w>0. 



(22) 



The ratio of the coefficients of the second derivatives is the propagation velocity squared of 
the spin-2 mode. Therefore, in terms of physical time, t = J a{T)dT and Planck units 7 = 1/2, 
the condition that a background have a real and non-superluminal spin-2 propagation speed is 



0<cl 



1 + 8/ 



• <1> 



(23) 



l + 8Hf 

where a dot indicates a derivative with respect to t, and H = a/a. Further, the condition that 
the spin-2 mode not be a ghost is 



< 1 + 8Hf, 



(24) 



The same strictures apply to the scalar (spin-0) mode, whose kinetic term is much more 
difficult to evaluate. Fortunately, our Lagrangian is a special case of a class of theories studied 
in Ref. [64]. There it was shown that the gauge-invariant combination $ = y? — H6(f)/(f) (one of 
the Bardeen scalars) satisfies 



where 



M («'«*) -4*-' 



Q 



96/72/2(1 + sHf) 



(1 + 12Hfy 
32H^P 



[1 + UHff 



-8f + 8Hf + {l + 8Hf)\3 + A^^^ 



H 



(25) 



(26) 
(27) 



If the spin-2 propagator is well behaved, the sign of the scalar time derivatives is also correct 
thanks to Eq. (24), and there are no spin-0 ghosts. Finally, the condition that the scalar 
propagation speed be real and non-superluminal reads 

0<clJl^'-4,]-'-^l^^<l. (28) 

- ° \^ sny 3(i + 8Hf) - ^ ' 

Having established these results, let us now focus on the particular class of models in hand. 
We have 

^ " 24 '^ ' ^^^^ 

where A = 24n ^^ /i^"+^. Without loss of generality, we fix to unity the coefficient of the square 
of the Riemann tensor in the lagrangian, i.e. 03 = 1, and note that the equation of motion, 
in a ffat FRW background, then yields 

(j) = Rl^ = 24H^- . (30) 

Since we are interested in the behaviour of these actions in the universe at late times, we study 
the attractor solutions for the CDDETT model when the matter is diluted away because of the 
expansion. In general [26] the attractor solutions can be written as a{t) oc V\ with p > 1. For 
the Gauss-Bonnet combination, it was found that the relevant accelerating power-law attractor 
is given by p = 4n + 1 . We see that 

with > for an accelerating universe, and therefore 

/ = 4(n + l)(p-l)^^, (32) 

A rfi An 4- 3 
/ = 96(n+l)(4n + 3)(p-l)^^^ = ^/. (33) 

Thus, it is clear that, if A > or ^^ = 1, (22) are satisfied for an accelerating universe (p > 1). 
However, the spin-2 propagation speed is 

_ l+8(4n + 3)/A 
1 + Spf/t 

which for the pure GB theory approaches 1 from above, and so the graviton propagates faster 
than light. 



For the scalar, condition (28) amounts to 



2 ,_, 4 V 8(p-4n-3)//t 



0<Cci= [l + -— r^^^^ < 1 . (35) 

- ° V 3^; 3 i+8pf/t - ^ ^ 

A sufficient condition to satisfy these relations is 

max(4, n + 7/4) < p < 4n + 3 , (36) 

which is clearly satisfied for the pure GB theory for any n > 1. 

To summarise: there are no ghosts or instabilities for the Modified GB attractor solutions 
for any n > 1[67]. However, the graviton propagates superluminally, which may render the 
pure GB theory inconsistent. 

Finally, the pure GB combination is not phenomenologically viable. In this case, (p oc H^'d 
and vanishes as a approaches to zero. This means that, for this combination, in the CDDETT 
model a = is a singularity of the equations of motion, and it not possible to change the sign 
of a. The Universe can never change from deceleration to acceleration. (This is reminiscent of 
the problem encountered in some other modified gravity theories [65].) 

B. Modified Gravity with matter 

In this section we study the propagation of the spin-2 and spin-0 modes when matter is 
present, which changes the background cosmology and hence the coefficients of the second 
derivatives in the action. 

As we discovered in the previous subsection, a realistic model should have ai ^ a^ so that 
the modification is not pure Gauss-Bonnet. The conditions (23,24,23,28,28) can be generalized 
in a straightforward but lengthy way (see [66]), and depend on the parameter b = ai/a^ — 1 
(recall we are assuming that a2 = — 4a3 as required to cancel fourth-order derivatives). The 
correct signs in the spin-2 propagator are assured if 

1 + AbfR + 8/ > 0, (37) 

l + 4bfR + 8Hf > 0, (38) 



with the propagation speed condition reading 



l+4bfR + 8Hf 
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The coupling of /(0) to the Ricci scalar outside the GB combination gives an extra propagating 
scalar degree of freedom \I' (see [66]), the other Bardeen scalar. The two scalar modes have the 
same speed of propagation Cg and, once again, if the graviton propagator is well behaved, there 
is only the following extra condition to be satisfied 

where 

Qi = 4b{fR + fR) + 8fH^ and Q2 = I + 4bfR + 8Hf . (41) 

These conditions put additional bounds on the parameter space spanned by the /x, Oj which 
define the CDDETT Modified Gravity theories, if we require that they hold at all times during 
the evolution of the universe. In particular, we may require that they hold when the Universe 
expands with a power-law, a oc t^, for which 

= ^ (42) 

R=^ (43) 

/ = r3t^"+^ (44) 

where 

Ti = 12p'^[3b{2p-lf + 2p{p-l)] (45) 

Ts = 6p{2p - 1) (46) 

n /7''"+2 



Therefore we have that 



2 l + 4r3t4"+2[6r2 + 8(n + l)(4n + 3)] 



^2 i + 4r3t4"+2[5r2 + 8p(n + i)] ' ^^^^ 

2 ^ ^ _ 32 p{n + l) 8 (n+l)(4n + 3-p) 

^° 3 6r2(4n + 2) + 8j92(n+l) 3 [612 + 8j9(n + 1) + l/4r3t4«+2] 

This relation holds for any power-law behaviour, unless either Fi or r2 becomes zero. We 
expect power-law expansion at early time, where the Universe should behave as an ordinary 
Friedmann model, and at late time where it approaches an accelerating attractor solution. 

Let us consider a matter-dominated Universe, p = 2/3, assuming a positive F3 to ensure the 
no-ghost conditions are satisfied. In this case it is clear that the tensor modes are superluminal, 
and their speed tends to one as t ^ 0. 

11 



On the other hand, for a late-time attractor solution, the accelerating power law is given by 



2(n + l)(4n + 1) - 3a + J'^n^a^ - 4(n + l)(4n + l)(5n + \)a + 4(4n + If in + 1)2 ^ ^ 

^ = w^) ^^^ 

where a = (12ai + 4a2 + Aa^) / {12ai + 3a2 + 203) = 66/(66 + 1). The exponent p is real, for 
large n, if a < | n or a > 8 n. Furthermore at late times, Eq. (48) takes the approximate form 



^2_ &r2 + 8(n + l)(4n + 3) 
^2 "^ 612 + 8p(n + 1) 



(51) 



which means that for non-superluminal behaviour we require 

p > 4n + 3 . (52) 

It should be noted that the value a = 1 cannot be considered, because this would imply that 
02 7^ — 4:0,3- For large n, one requires 

to avoid superluminal tensor modes. 

At late times, Eq. (49) takes the approximate form 

2 32 p{n + l) 8(n + l)(4n + 3-p) 



" 3 6r2(4n + 2) + 8p2(n + l) 3 feEs + 8p(n + 1) 

Eor large n, Eq. (54), implies that the scalar modes are not superluminal if 



(54) 



24 
-y<«<l- (55) 

For large n, in the region —24/5 ^ a ^ —8/5 both tensor and scalar modes are not superlumi- 
nal. In the case of n = 1 one can see that the allowed region is —3.793 ^ a < — ^ ~ —1.965, 
where the lower bound is given by c^ = and the upper one by C2 = 1. This additional 
constraint would rule out the region 0.9 < a < 1 (where p is real) identifed by Mena at al. as 
producing cosmologies consistent with Supernovae data. 

A plot of the allowed values of a and n, after imposing the no-ghost constraint, or the 
no- superluminal constraint or both, is shown in Fig. IIIB. 

IV. CONCLUSIONS 

The search for a satisfactory model for the acceleration of the universe has been pursued 
in many different ways. Recently, models attempting to explain such behavior by changing 

12 



n 




FIG. 1: Contour plot in the {a,n) plane for the constraints. The light grey area corresponds to the 
region in which only the no-ghost constraint (Eq. 38) holds. The darker area represents the points 
at which both the no-ghost and the positive-squared-velocity conditions (Eqs. 37, 38, 40) hold at the 
same time. Finally the darkest region is the region of the plane at which all the constraints (no-ghost, 
< c^ < 1, Eqs. 37, 38, 40, 39) hold for both scalar and tensor modes. 



the gravity sector have been proposed [25-37]. In particular, the CDDETT model [26], has 
the attractive feature of the existence of accelerating late-time power-law attractors, while 
satisfying solar system constraints [53]. 

In this paper we have investigated the consistency of the propagating modes (tensor and 
scalar) for the action 



S 



d X 



-9 



iR + e, 



/i 



4n+2 



(56) 



(ai i?2 + a2 Ra(3 R'"^ + as Ra^yS R""^^^)' 

In order for this action to be ghost-free, it is necessary but not sufficient, to set 02 = —4 03 
[39, 53] so that there are no fourth derivatives in the linearised field equations. What remained 
was the possibility that the second derivatives might have the wrong signs, and also might 
allow superluminal propagation at some time in a particular cosmological background. For 
example, for the case Oi = 03, for which the modification is a function of the Gauss-Bonnet 
term, we found that the accelerating power-law attractor solutions give propagators with the 
correct signs, but with a spin-2 mode propagating faster than light. 

We have also examined the general second order CDDETT Modified Gravity theory in a 
FRW background with matter, which is parametrized by the energy scale /i and by 6 = 01/03 — ! 
- the deviation of the Ricci scalar-squared term from that appearing in the Gauss-Bonnet 
combination - or equivalently a = 66/(66 -|- 1). We found that the theories are ghost-free. 



13 



but contain superluminally propagating scalar or tensor modes over a wide range of parameter 
space. In conclusion, we note that there are likely to be further constraints from compatibility 
with CMB data as we have changed gravity on large scales quite significantly. To investigate 
this point is beyond the remit of the current paper. 
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APPENDIX A: SCALAR FIELDS AND MODIFIED GRAVITY 

In this Appendix we demonstrate how actions for non-standard models of gravity can be 
rewritten in the form of Einstein gravity with a non-minimal coupling to one or more scalar 
fields. This is not a new result (see most recently e.g. [39]), but the Modified Gravity action 
(2) is a special case which needs separate consideration. 

Consider the action 

S= fd^x./^F{Z,) , (Al) 

where Zi are monomials in the curvature invariants, with Zi, Z2 and Z3 defined earlier, but 
where we allow the possibility of higher order terms. If the function FiZi) is at least twice 
differentiable (except possibly at isolated points), then this is easily seen to be equivalent to 
the action 

S4,. = jd'xV^ [F(0) + {Z, - 0.)F,(0)] , (A2) 

where 0j are a set of auxiliary scalar fields, one for each of the terms Z^, and Fi = diF. The 
first variation is 

5% = J d^x v/^ [{Z, - <Pi)F,,{<P)6<P, + (5Z,F,(0) + ]^C^,5gl (A3) 

where C^^ = F{(j)) + {Zi — (f)i)Fi{(f)). We immediately see that, provided the matrix of second 
derivatives Fij = didjF{(f)) is non-singular, 0j = Zi, and we return to the original action (Al). 
This was one of the results in [39]. 

However, the possibility of a singular matrix was not considered. In models of the form (2), 
there exist degeneracies in the parameters of the form 

F{Z, - AX) = F{Z,) , (A4) 

where A^ are arbitrary constants, and f ° are orthonormal vectors in the space of curvature 
invariants. In this case, ff-Fj = 0, -F^f ? = 0, and higher derivatives also vanish. 
The solution of 6S/6(f)i = is now 

(Z, - 0,) = A^ , (A5) 
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with Aa again arbitrary constants. By substitution we find 

Sa = J d'x v^ [F{Z, - A„<) + A,<9,F(0)] = J d'x V^F{Z,) , (A6) 

which is again our original action. 

As our Modified Gravity example suggests, we can reduce the number of scalar fields by 
taking linear combinations of the Zi normal to the subspace spanned by the f f . If this subspace 
is spanned by wf, then we can define a new set of variables Z'j^ = wfZi, Z'^ = vfZi. The 
function F is now independent of Z^, so we can write 



S[g,u, U = j d^x^/^g [^(0:4) + {Z'^ - <l>'A)FAm . (A7) 

In the CDDETT model, the Lagrangian density depends on Zi, Z2, and Z^ only through the 
combinations Zi, —4^2 + Z^. There is one degeneracy, and so there are only two scalar fields 
required to put the action into the linearised form of Eq. (A2). When 6 = in (6), the field 
associated with Zi = R^ may be trivially solved for to give the Einstein-Hilbert term. 



